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Abstract. A variety of physical phenomena can lead to viscous effects. In this 
talk we review several sources of shear and bulk viscosity with an emphasis on 
the bulk viscosity associated with chiral restoration. We show that in the limit 
of a second order phase transition, the viscosity peaks in a singularity at the 
critical point. 
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1. Introduction and Theory 

Viscosity has attracted remarkable attention at RHIC during the first years of run- 
ning. In particular, experimental observations of large elliptic flow have pointed 
to a small shear viscosity and inspired the term "perfect liquid" [^. In this talk, 
we review the general theoretical definition of viscosity, then show how five differ- 
ent physical effects can lead to non-zero viscous coefficients (Sec. 2). We focus 
on bulk viscosity, and show that one can find large, even singular, effects in the 
neighborhood of T^ (Sec. 3). Although the present study focuses on understand- 
ing the behavior and physical explanation of the coefficients, we speculate on the 
experimental manifestations large viscosities might bring about. 

In non- viscous hydrodynamics the elements of the stress-energy tensor depend 
only on the energy density e and particle-number densities n when viewed in the 
rest frame of the matter. 



T, 



non. Vise 



(r,t) = <5yP(e(r,f),r?(r,t)), 



(1) 



where e and n are implicitly functions of r and t. The tilde denotes that Ty 
is evaluated in a frame where the collective velocity u(r) = 0. In Navier-Stokes 
hydrodynamics, viscosity is incorporated by altering T so that it includes terms 
proportional to the velocity gradients dui/drj. 



T, 



(N.S.) 



r 7^/ ^N / ^N i' dui duj 2^ 



u6i 



S(e,n)(V-u)<5y, (2) 
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Here, 77 and B are the shear and bulk viscosities. Since V • u = {de/dt)/{P + e), 
the bulk viscosity can be interpreted as a correction to the pressure due to the 
changing energy density, whereas the shear viscosity describes the asymmetry of 
Tij due to an anisotropic expansion. In non-viscous hydrodynamics accelerations 
are proportional to the gradient of the pressure, while in general, accelerations arise 
from derivatives of the stress energy tensor. 



, duj d ^ 

dt dxi 



(e%+^«.)^--7T-T'.,- (3) 



Thus, the components of the stress-energy tensor can be considered as representa- 
tives of the pressure in a given direction, and any reduction/rise of Tij from viscous 
effects will result in a slowing/acceleration of the expansion in that direction. 

Viscous coefficients can be expressed in terms of correlations in the stress-energy 
tensor through Kubo relations. These are derived by considering alterations of Tij 
due to a perturbation V in linear response theory, 

{dT,,{r = 0)) = -{t/h) I d\'m,{r = 0), V{r')]), V{r') = r'Mu,)To,{r'), 

(4) 
where the perturbation represents the change to the Hamiltonian due to boosting 
according to a linear velocity gradient. After integrating by parts and applying the 
conservation of the stress energy tensor, dtTi^ = —djTij, one can derive the Kubo 
relation, 

6{T,j{r - 0)) - {i/h) f dV'r(,([Ary(r = 0),ATki(r')])dkUi. (5) 



Here, AT^i refers to the difference with respect to the average of Tki a,t t = —00. 
For i y^ j, symmetries constrain k and I to equal i and j, which allows one to read 
off the shear viscosity from Eq. (0) , 

ri = {^/h) f d\'r',{[AT,,{0),AT,,{r')]), ^^J (6) 

Jr'„<0 

= i^-So 2^ / d*r'e^-''{[m, (0), AT,, (r')]). 

By considering the case where diUj — (l/3)(5y V • u, one can inspect Ta in Eq. ((SJ 
to find the bulk viscosity, 

B = mh)Y,f dVr^([AT,,(0),Ar,,(r')]) (7) 

, Jr'„<0 
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The Kubo relations, Eq.s © and iQ, are fully consistent with quantum mechanics. 
The classical limit can be obtained by applying the identity [|2j, 

e-^"Vit)^e'^^^W{t)e-f''', (8) 

to one of the terms in the commutator in Eq.s © or (fTJ, then keeping the lowest 
term in h, 

Tr e-''^[Ar,,(0), ATki{r)] « -thpTi dt {h.T,,{0)DeltaTki{r)) , (9) 

which after an integration by parts gives the classical limit of the Kubo relations, 

r^ ^ P [ dV(AT,,(0)AT,,(r')), t^j (10) 

Jr',<0 

B « {vmJ2[ d\'{ATuiO)AT,,{r')} (11) 

Although the Kubo relations are difficult to interpret physically, they do make it 
clear that viscosity is related to the size and to the damping of fluctuations of the 
elements Tij. If fluctuations in T^ (at fixed energy) are large, or if they are slow to 
relax, a large viscosity will ensue. 

2. Five Sources of Viscosity 

Viscous effects arise whenever the elements of the stress-energy tensor, Ty, have 
difficulty maintaining the equilibrium values due to a dynamically changing system, 
i.e., one with velocity gradients. In this section we briefly review five physical sources 
of viscosity, the first three of which have already been explained in the literature. 

1. Viscosity from non-zero mean free paths: This is the most commonly 
understood source of viscosity. It is straight-forward to see how a non-zero 
collision time leads to an anisotropy for Tij by considering a velocity gradient 
for a Bjorken expansion, Uz — z/t, i.e., dzUz — I/t^OxU^ = dyUy = 0. We 
consider a particle of momentum whose momentum is p^ (t) when measured 
in the frame moving with the collective velocity corresponding to its position. 
In the absence of collisions, p'z will fall with r since the particle will asymp- 
totically approach a region where its velocity equals the collective velocity, 
p'z{t + St) = p'^{t)t/{t -f 5t). Meanwhile, p',j, and p' are frozen. The result- 
ing anisotropy in the stress-energy tensor is easy to derive, which yields the 
following expression for the shear viscosity [|3j, 

7? - (4/3)Pt„ (12) 

where Tc is the collision time. It is also easy to see how such an expansion does 
not yield a bulk viscosity for either ultra-relativistic or non-relativistic gases. 
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In those cases an isotropic expansion scales all three momenta proportional 
to l/r which maintains thermal equilibrium, and collisions do not play a role. 
This is not the case when m ^ T, or especially if the gas has a mixture of 
relativistic and non-relativistic particles. 

2. Viscosity from non-zero interaction range: If the range of interaction 
between two-particles extends a distance R, interactions will share energy be- 
tween particles from regions with different collective energies. A particle at 
r = 0, where the collective energy is zero, will share energy with particles 
whose collective energy is {l/2)m{Rdru)^. For Boltzmann calculations, the 
viscosity will be proportional to PR'^/tc [0], with the constant of proportion- 
ality depending on the scattering kernel. Both bulk and shear terms result 
from non-zero interaction range. In Boltzmann calculations, the range of the 
interaction can approach zero for fixed scattering rates if the over-sampling 
ratio is allowed to approach infinity. Although this solves causality problems [ 
|S] , it simultaneously eliminates viscous terms arising from finite-range scatter- 
ing kernels, which might or might not be desirable. This has profound effects 
on calculations of elliptic flow, which can vary by a factor of 2 depending on 
the range of the scattering kernel [ ::4 . 

3. Classical Electric Fields: Color fiux tubes form after the exchange of soft 
giuons between nuclcons passing at high energy, and might also be formed 
during rapid hadronization. Additionally, longitudinal color electric fields 
might be created during the pre-thermalized stage of the collision (color-glass 
condensate). Since these fields tend to align with the velocity gradient, they 
can be a natural source of shear viscosity. In fact, if the fields are purely 
longitudinal, the elements of the stress-energy tensor become T^z = — e, T^x = 
Tyy = e. Thus, the transverse pressure becomes three times as stiff as a 
massless gas, P = e/3, which is usually considered a stiff equation of state. 
The negative longitudinal pressure signifies that the energy within a given 
unit of rapidity is increasing during the expansion, similar to the negative 
work associated with stretching a rubber band. 

4. Non-equilibrium chemistry: Chemical abundances can not keep up with 
the expansion if the rate at which equilibrium abundances change is not much 
smaller than the chemical equilibration rate, l/rchom, 

^ = -(l/rchc„0(^-^cq). (13) 

If the equilbrium number is slowly changing, the abundances will vary from 
equilbrium by an amount, 

5iV = -r,hc„.^. (14) 
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To associate this departure from equilibrium as a viscosity, one must consider 
the corresponding change in the pressure, 



. _ dP 
dn 



fixed . dt ds 



where the second relation exploits the fact that entropy is conserved in a slow 
expansion. The bulk viscosity is then. 



_ dP 

dn 



fixed e ^* 



'^. (16) 



The bulk viscosity will be large whenever the equilibrium number is rapidly 
changing, e.g., the temperatures are falling below the masses, or masses are 
rising due to restoring chiral symmetry. If the hydrodynamic equations ex- 
plicitly treat particle numbers as current obeying chemical evolution rates, 
chemical non-equilibration would not need to be accounted for through vis- 
cous terms. 

5. Viscosity from dynamic mean fields: Bosonic mean fields, such as the 
a field, obey the Klein-Gordon equation. For fluctuations of wave number 

^Aa(i) = -(m„(T)2 + fc2)Aa(t)-r|Aa(t), (17) 

A<T(i) = o-(i)-acq(e), 

where <Toq(e) is the equilibrium value of the condensate which is non-zero for 
fc = 0. The value of a is determined by minimizing the free energy, while the 
mass is related to the curvature of the free energy near the minimum, 

^F(a,r)=0, ml^^{T)^^F{a,^,T). (18) 

One can see the equivalence of Eq. (|17|) with the differential equation for the 
harmonic oscillator after performing the following substitutions, 

fch.o./wh.o. ^ m^, 7h.o./™h.o. ^ r, (19) 

where 7h.o. is the drag coefficient for the harmonic oscillator, fch.o. is the spring 
constant and Wh.o. is the particle mass. For the harmonic oscillator, the mean 
value of the position x is altered if the equilibrium position is moving. The 
amount of the change was consistent with the drag force "fdxoq/dt being equal 
and opposite to the restoring force kSx. The corresponding result can be 
derived for the damped Klein-Gordon equation, 

c 7h.o. dXgq i "O'cq /r)n^ 
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where 6a is the mean offset from the equilibrium value. Thus, rri^ determines 
the restoring force, while F describes the drag. Finite-size effects could be 
estimated by replacing vi^ with ni? + fc^, where k^ would be given by the 
finite size, k ~ 1/ L. The resulting bulk viscosity is. 



B = 



dP 



F da, 



cq , 



fixed 



TO^ 



ds 



(21) 



The bulk viscosity is then large for energy densities where a is rapidly varying, 
or for when nif, is small, i.e., the critical region. 



3. Bulk Viscosity in the Linear Sigma Model 

Both of the last two sources of viscosity described in the previous section can be 
of special importance during the chiral transition. First, since masses are changing 
suddenly near Tc, chemical abundances should easily stray from equilibrium. Sec- 
ondly, the mean field, which is zero above Tc suddenly changes, and given the small 
masses in this region large bulk viscosities are expected. 
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Fig. 1. For the linear sigma model, the sigma field and mass are shown as a 
function of the temperature in the left panels. Near Tc, the masses fall to zero and 
the mean value of the field changes rapidly, which gives rise to a sharp peak in the 
bulk viscosity. The pressure and temperature are displayed in the right-side panels 
as a function of energy density. 



As an example, we consider a simple example of a linear sigma model, where 
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Fig. 2. For a Bjorken expansion (V • u = 1/t), the pressure is plotted alongside 
Tii — P — BV ■ u, to demonstrate the significance of viscous terms near Tc. Viscous 
terms are larger and sharper for couplings close to the critical coupling. 



the coupling of the sigma model to the quarks provides the quark mass [El El, 



iJ = — ctW 
2 



\4 

— {a^ - fl+ ml/\^f - hga + i/quarks(m = ga), (22) 



assuming only up and down flavored quarks. The resulting equation of state and 
values for ma and a are displayed in Fig. ^for A^ — 40. For couplings g < gc = 
3.554, the transition is a smooth cross-over, while for g = gc the transition is second 
order, and for g > gc & first-order phase transition ensues with Tc = 172 MeV. From 
Fig. n one can see that rria becomes small in the same region that the field rapidly 
changes, thus one expects a peak in the bulk viscosity. 

The bulk viscosity was calculated according to the methods described in (4) 
and (5) of the previous section assuming that the width F = 400 MeV, and that 
the chemical equilibration time Tchom=l fm/c for a density of one quark per fm'^ 
and scaling inversely with the density. For a Bjorken expansion V • u = l/r, and 
assuming an isentropic expansion starting with e = 8 GeV/fm'^ at r = 1 fm/c, we 
calculated both P and B as a function of r. To illustrate the size of the effect, we 
display both P and the Navier-Stokes expression Tu = P — B\7 • u as a function 
of the energy density for a Bjorken expansion in Fig. [3 The effect is certainly 
dramatic, especially for g « gc- However, the Navier-Stokes expression is only 
applicable for small expansion rates. We expect Israel-Stewart [|H1 IHI equations for 
hydrodynamics to result in moderated effects compared to Navier-Stokes, though 
they should give identical results for small velocity gradients. 



4. Summary 

The simplicity of the Kubo relations, Eq.s (jBJ and Eq. ((Tl, masks the wide variety of 
physical sources of viscosity. The one common aspect of the various sources is that 
non-zero equilibration times, non-zero interaction ranges can always be identified. 
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In this talk, we focused on bulk viscosities associated with the chiral transition. In 
general, one would expect such effects whenever a system needs to rapidly rearrange 
its basic structure. In this sense these effects have much in common with super- 
cooling or hysteresis. 

The implications for dynamics should be that the matter accelerates more 
quickly due to the higher gradients in Tx^x when the interior energy density is above 
the critical region. However, once the matter flows into the viscous region of energy 
densities, there should be a slowing down and a reduction of surface emission. This 
trend would be in the right direction to explain HBT measurements which show a 
rapid expansion with a sudden disintegration [I10|. but the potential magnitude of 
the effects are not yet known. 

Finally, we re-emphasize that if one were to solve for the evolution of the mean 
fields or chemistry alongside solving the hydrodynamic evolution equations, one 
might be able to neglect some of these effects. If these effects are large, the proper 
conclusion may be that rather than absorbing these effects into viscous hydrody- 
namics, one should treat non-equilibrated degrees of freedom more explicitly. 
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